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We determine the charge radius of the proton by analyzing the published low momentum transfer electron- 
proton scattering data from Mainz. We note that polynomial expansions of the form factor converge for mo¬ 
mentum transfers squared helow 4m^, where myc is the pion mass. Expansions with enough terms to fit the data, 
hut few enough not to overfit, yield proton radii smaller than the CODATA or Mainz values and in accord with 
the muonic atom results. We also comment on analyses using a wider range of data, and overall obtain a proton 
radius Re = 0.840(16) fm. 


I. INTRODUCTION 


Much remains to be learned about the proton. After a half- 
century of study, we still do not know what its size is, where 
its spin comes from, and how its mass is generated from light 
quarks and gluons. Particularly troubling is the matter of 
the proton’s charge radius. Re- This was first measured to 
be approximately 0.8 fm by Hofstadter and collaborators in 
the 1950s via elastic electron scattering m. The value of 
Re has been steadily refined over the years through electron 
scattering and hydrogen energy level measurements, recently 
reviewed in Refs. ma. The CODATA group ilia, using 
available electron-based data through 2014, quotes a com¬ 
bined value of Re — 0.8751 ±0.0061. The recent electron 
scattering experiment in Mainz EHH, which quotes a value of 
Re = 0.879 ± 0.008, is included in this CODATA value. For 
many years Re had remained relatively stable at about 0.88 
fm, until this value was called into question by Lamb shift 
measurements in muonic hydrogen, which yielded a value of 
Re = 0.84087 ± 0.00039 fm m [H. This radius is 7 stan¬ 
dard deviations away from the CODATA value. The proton 
size puzzle leaves us with three options: the hydrogen Lamb 
shift and elastic electron scattering experiments have erred in 
extracting Re, the muon Lamb shift measurement is precise, 
but inaccurate, or there is new physics that affects the muon 
differently than the electron, rendering the theory behind the 
muonic Lamb shift calculations incomplete miHisi- 

In this paper, we explore whether the published, high- 
quality ep elastic scattering data at low-Q^ from Mainz could 
be consistent with the muonic Lamb shift determination of 
Re- Extracting Re from elastic electron scattering is as sim¬ 
ple, or as difficult, as measuring the slope of the electric form 
factor Ge{Q^) as a function of the squared four-momentum 
transfer as goes to zero. However, since the differen¬ 
tial cross section diverges at small scattering angles (low Q^), 
these measurements are very sensitive to beam alignment and 
angle determination. No ep measurement extends to = 0, 
although some get close. Mainz currently holds the record, 
with measurements at as low as 0.0038 GeV^. This mod¬ 
ern data set, with 1422 data points in the range from the lower 
limit to about 0.98 GeV^, is the best, most precise, and most 
extensive available. Therefore, it is the Mainz data that we 
choose to explore. 

The charge radius is given by the second term in the expan¬ 


sion of the electric form factor, 

Ge{Q^) = 1-\rIQ^ + c2Q^ + ---- ( 1 ) 

o 

Using data at very low Q^, one can hope the term, the 
curvature term, is small so that the charge radius Re can be 
determined without having to model the shape of Ge over a 
wider range of In the past, the size of the uncertainties on 
the data at very low has meant that one could not extract 
an accurate charge radius without extending the data range to 
include not-so-low The term then can become notice¬ 
able, depending on how far the data range is extended. There 
was an early example of Simon et al. data iflTl where the fitted 
coefficient of the Q‘^ term was small, albeit with large uncer¬ 
tainty, and the extracted charge radius was also small. Us¬ 
ing this example, some workers {e.g., ifTSl ) advocated includ¬ 
ing data at still higher ^ to obtain a larger curvature. This 
also led to a larger extracted proton radius. Including higher 
2^ does not have to mean including data at all available Q^, 
and we have the example of Ref. m using only data with 
< 0.62 GeV^ (2 < 4 fm-'). 

Sick and Trautmann |j20l, in fact, suggest that data from 
0.014 to 0.056 GeV^ in (or Q from 0.6 to 1.2 fm-') is 
most crucial for finding the proton radius. The consideration 
of what data range is sensitive to the proton size, given the 
foregoing discussion, must depend on the accuracy and preci¬ 
sion of the data. As the data improve, the range of needed 
to obtain the proton radius will decrease. With the new Mainz 
data now available, we will explore the possibility that we can 
obtain a good proton radius result using only data with a rather 
low maximum Q^. 

The Mainz data HSU data enjoy state-of-the-art radiative 
and Coulomb corrections. One of the three spectrometers 
was used as a luminosity monitor to control systematic un¬ 
certainties, and the other two spectrometers measured sepa¬ 
rate kinematic points simultaneously. The data are dominated 
by point-to-point systematic uncertainties from background 
subtractions, drift-chamber inefficiencies, normalization fac¬ 
tors, angle determinations, and the afore-mentioned correc¬ 
tions. The slight leeway we have in fitting these data is to 
make small rescalings of the 34 normalization sets in the ex¬ 
periment and to enlarge the point-to-point error bars if the 
fluctuations of the data indicate that the quoted uncertainties 
are too small. We limit our form factor fits to the Mainz data 
set because the systematic differences between separate ex¬ 
periments can introduce systematic effects in global fits, and 


2 


the 1422 Mainz points already dominate the sample of world 
data below = 1 GeV^. 

We shall advocate for fits using the 243 Mainz data points 
with below 0.02 GeV^. In addition to the general principle 
that using only very low will free us from model depen¬ 
dence incurred in extending the fits to higher Q^, we are also 
using only data from a limited number of spectrometer set¬ 
tings. This substantially frees the data from any drift arising 
from the normalization adjustments that reconciled data from 
different spectrometer settings, and averts overfitting of in¬ 
flections or statistical fluctuations in the data that can happen 
when using fit functions that contain many fit constants. 

In the following, Sec. [represents the formalism pertinent to 
our discussion; Sec. [represents the extractio n an d discussion 


of the proton radius using low data; Sec. IV presents a fit 
based on the full range of the Mainz data; Sec. ^highlights 
our final results; and Sec. VI gives our conclusions. 


II. FORMALISM 


The dipole form factor 
Gd = 


(l+GVO -71 GeV^) 


2t2’ 


( 8 ) 


suitably normalized, has been used for many years as a bench¬ 
mark approximation for both Ge and Gm- The Mainz group 
presents data on the cross section a as its ratio to the cross 
section calculated using the dipole form factors. Go, 


a 

Od 


eGl -f xGl 
eGl + xii^Gl 


( 9 ) 


From this. 


Ge{Q^) =Gd{Q^) 


(- 


1/2 


l + TjJ. 


2Gm/(MpG£)2- 1 


e + Tfip 


-1/2 


( 10 ) 


When an electron of energy E scatters from a proton at rest 
through an angle 0 and exits with energy E', the 4-momentum 
transfer squared is, 

= -q^ =AEE'fim^^. (2) 

In the Born approximation, the ep elastic scattering cross 
section can be written in term of the electric and magnetic 
Sachs form factors, Ge{Q^) and Gm{Q^), 


da 

/da'' 


dQ. 

\dQ., 

^Mott (1 + '^) 


The Mott cross section is 


da\ 


gUq^) + -gUq^) 


4a"^cos^ I E'^ 


Mott Q 

a is the fine structure constant. 


E ’ 


e = ^1 4-2(1-f T)tan^ , and 

Q2 4^2 ’ 


( 3 ) 


( 4 ) 


( 5 ) 


where v is the energy transferred by the virtual photon and M 
is the proton mass. Further, 


1-f (2£//W)sin2f ■ 

The electric and magnetic form factors at = 0 are nor¬ 
malized to correspond to the nucleon charge in units of e and 
the nucleon magnetic moment in units of the proton magneton 
Pn = el{2M), such that 


G£(0) = 1, and Gm(0)=Mp« 2.793. (7) 


For < 0.02 GeV^, the quantity in square brackets above, 
for reasonable values of the Ge /Gm ratio, differs from unity 
by no more than 140 parts per million, and plays no significant 
role in the extraction of Ge- For general Q^, we most often 
will obtain Ge using the Ge /Gm ratio obtained from recoil 
polarization experiments, mostly at higher Q^. The data from 
JLab indicate that at least for < 8 GeV^, 


Ge ^ 

~ 8 GeV2 ■ 


( 11 ) 


Returning to low one can use the radius expansion. 


Pp^ = l-URl-Ri)Q\ ( 12 ) 

(JM D 

where Rm is the magnetic radius and 


Ge 



1-f 


nlQ^Rl-Rj 

12M^ e + Tpj 



( 13 ) 


Specifically, Bernauer et al. Ih) obtain Re = 0.879 ± 0.008 
fm and Rm — 0.777 ±0.017 fm from their fits. The latter is 
significantly smaller than most other fits, where typically Rm 
is similar to Re, but none-the-less, the extraction of Ge from 
the data at low is unaffected by the spread of suggested 
values for Rm even at the high level of accuracy needed for 
the present investigation. 


III. ANALYSIS FOR < 0.02 GEV^ 


Looking only at data for below 0.02 GeV^, there are a 
plethora of data points from the Mainz experiment. In this 


range, the term in brackets in Eq. 10 is, for all practical pur¬ 
poses, unity. Mainz has 243 data points for < 0.02 GeV^. 
Limiting consideration to only spectrometer B to minimize 
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FIG. 1. (Color online) Linear plus quadratic fit to all Mainz (2010) 
data with < 0.02 GeV^. 


cross calibration uncertainties, still leaves 209 data points. Al¬ 
though the overall normalization may be uncertain by 1 - 2 %, 
any relative systematic uncertainties that could lead to a false 
dependence are thought to be small. 

We can fit the data using a linear plus quadratic in form 
for Ge, 


GEiQ^)=ciil+C2Q^ + C3Q'^), (14) 

where C 2 = —/?|/ 6 . Using all available points below 0.02 
GeV^ gives the result shown in Fig. The per degree of 
freedom (dof) for the fit is 1 . 00 , the normalization constant 
is 0.9992 ±0.0003, Re = 0.850 ±0.019 fm, consistent with 
the Lamb shift results, and C 3 = 4.5 ±5.6 GeV^^. The central 
value of C 3 is positive as one might expect from a nonrelativis- 
tic expansion of the form factor, but is statistically consistent 
with zero. 

Without prior expectations, one could use a rule of thumb 
for fitting, namely to discard terms in the fit that do not im¬ 
prove the x^/dof. A linear fit, Ge{Q^) = ci(l +C 2 Q'^), leads 
to the same j^/dof = 1.00, and values ci = 0.9986 ±0.0003 
and Re — 0.835 ± 0.003 fm (using the diagonal term in the 
error matrix). From a statistical viewpoint, a radius as large as 
0.88 fm is unfavored. Flowever, the correlations are such that 
the more positive the curvature, the larger the extracted proton 
radius. 

Fits like the ones just presented have been criticized be¬ 
cause of expectations that the curvature could be larger than 
the apparent results from the quadratic fit and certainly not 
zero as in the linear fit. To investigate the effects of curvature 
when fitting a low-Q^ data set, we expand the electric form 
factor as 


GEiQ^) = l-lRlQ^ + ^RtQ^-^R^Q^ ( 15 ) 

where the coefficients are suggested by nonrelativistic mod¬ 
els, where Re is the rms proton radius, = (r^), and 
b 2 = (r'*)/(r^)^, and b^, = (r®)/(r^)^. The coefficients can be 
calculated using exponential, Gaussian, and uniform model 


charge distributions p(r) = poe p(r) = poe and 

p(r) = po0(c — r). For these three cases, respectively, b 2 = 
5/2, 5/3, and 25/21, and bs = 35/3, 35/9, and 125/81. 
The fits using < 0.02 GeV^ data yield Re = 0.859(3), 
0.851(3), and 0.846(4) fm, respectively, each with a j^/dof 
= 1.00. One of these results is almost neutrally between the 
muonic and electronic radius values, while the Gaussian and 
uniform distributions even with the pre-chosen curvature term 
give results commensurate with the muonic Lamb shift value. 

The fit just discussed is one example. One may inquire what 
results follow with different 2 ^,^ (always with 2 mar < 4 to^) 
and different orders of polynomial. There will be two criteria 
for an acceptable fit. One is that the is low enough, on 
the order of 1 per degree of freedom. The other is that the 
highest order term in the polynomial is not well determined, 
as judged by the uncertainty limit on its coefficient, with the 
previous terms well determined. This will imply that the fit 
has omitted no important term, and is good. Fig. shows 
X^/dof and Re from a number of linear, quadratic, and cubic 
(in Q^) fits, with the small numbers indicating for s^^h 
example in multiples of 0.01 GeV^. 
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FIG. 2. (Color online) ;t;^/dof vs. Re for polynomial fits to the data 
set, showing fits with first, second, and third order polynomials in 
using data with 0.01 < 2max < d.OS, GeV^. The small numbers near 
each data point gives Q^„ax in multiples of 0.01 GeV^. A good fit 
should have a good x'^ and a sufficient but not oversufficient number 
of parameters, as further discussed in the text. The quadratic fit “2” 
and the cubic fit “4” satisfy the criteria. 

The linear fits cannot sensibly satisfy the second criterion 
and still give a radius, but they are included to show what hap¬ 
pens when Q^xix increases and the number of terms in the fit 
does not. The diagonal elements of the error matrix are aston¬ 
ishingly small, but the fits are poor, as judged by x^- Three 
quadratic fits are shown, with QI^^x indicated on the Figure. 
The lowest one is our prime example, which satisfies all crite¬ 
ria. It has an acceptable x^ and a small contribution from the 
actual term in the polynomial. The next still has an accept¬ 
able x^, but the coefficient of the term is stringently deter¬ 
mined (circa 10 % uncertainty), leading to worries that a cubic 
expansion is not flexible enough for this and so should 

not be trusted. The last point in this series sees some rise 
in and with tightly determined coefficients again indicat- 


• • linear in 

^ ♦ quadratic in 

A cubic in 
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ing underfitting at this For the two cubic fits shown, the 
one with lower Smax reasonably satisfies the criteria, while the 
upper one has its coefficient tightly determined, again in¬ 
dicating underfitting. Quartic fits did not give suitable results 
with Q^cix limited by the nn threshold. Some values, 
although happening to give an Re that matches the muonic 
hydrogen results, leads to fits with several poorly determined 
coefficients, and for other 2max ^igi^ pattern of the coef¬ 
ficients does not match the alternation expected for smooth 
charge distributions and the Fourier transform formula. The 
overall conclusion from further examinations of fits to the the 
low 2^ range of the data, is that when the criteria for good 
fits are satisfied, the proton radius accords with the muonic 
hydrogen value. 


for Ge{Q^), in which C 2 = {Re l^- A truncated contin¬ 
ued fraction is a ratio of polynomials, and it resembles Fade 
approximates. The continued fraction could be dangerous, be¬ 
cause it can have singularities in the spacelike region when¬ 
ever one of the constants c, is negative. However, if singulari¬ 
ties do not occur within the fit range, the continued fraction is 
acceptable, and it allows a wide range of shapes covering sev¬ 
eral orders of magnitude with relatively few parameters. On 
the other hand, fit functions with few parameters are unable to 
capture small inflections in the data. Since there is no theo¬ 
retical restriction that forbids inflections, many people believe 
that they should exist. We have looked for inflections in the 
data, and we find no persuasive argument to include them in 
our parameterization. The low-Q^ expansion of the continued 
fraction is 


IV. THE FULL RANGE 

We have so far concentrated on using low-Q^ data to obtain 
the proton radius, but there is interest in considering the full 
data set. There are several topics to discuss. Can a smooth 
function, with relatively few parameters give an acceptable fit 
to the data, and what is the radius that follows from such a fit? 
What is the value and effect of fitting with more parameters? 
With more parameters one can fit systematic deviations or sta¬ 
tistical fluctuations from what is really smooth data, which 
can mar the overall fit and skew the extrapolation to the pro¬ 
ton radius. Also, with more parameters, there is a tendency 
for extensions outside the fit region to rapidly deviate from a 
properly smooth continuation of the data. 

Another item to consider is that the polarization method for 
obtaining the form factors has shown that Ge falls relative 
to Gm with increasing momentum transfer, approximately as 
R-vGeIGm = 1 - eV(8 GeV^) EBEl. We used this earlier 
when obtaining Ge from the cross section data. For low 
data the difference between Ge obtained using the polariza¬ 
tion results and using scaling, jJ-pGE = Gm, is minor. How¬ 
ever, the full range of the Mainz data gives unexpected sup¬ 
port for the polarization result. This is surprising, considering 
that it is a Rosenbluth experiment without hard two-photon 
cotTections, while all earlier Rosenbluth results gave scaling 
(IIpGe/Gm « 1). Also in the absence of two-photon correc¬ 
tions, a reduced cross section at fixed should be linear in 
e (see Eq. 0). Two-photon corrections change the slope in 
e, and may also give some and higher dependence, which 
will be sought in the data in the ensuing subsections. 


A. Full Range Fit 


For our analysis of the full data set, we have chosen a con¬ 
tinued fraction (CF) form 


m^) 


Cl 


1 + - 




1 + 




1 + 




(16) 


/ iowq ( 2 ^) =Ci[1-C2Q^+C2{c2+C3)Q'^ 

+ C2iic2+C3)^-C3C4)Q^ +■■■]■ (17) 



FIG. 3. (Color online) Four-parameter continued fraction fit to the 
full Mainz data set using PpGg /Gm = 1 — 2^/(8 GeV^). 


We extracted Ge vs. using Eq. (10 1 and fit all data to 
a 4-parameter continued fraction form. Adding a fifth pa¬ 
rameter did not improve the fit, so we limited ourselves to 
4 parameters. The j^/dof is 1.6, which is somewhat high, 
and we shall have more to say about this later. However, the 
data are well-fit on average in all regions of Q^. Erom this fit 
we obtain Re = 0.8389 ± 0.0004. The uncertainty on Re is 
small because a constraining fit form introduces information 
into the problem, in this case a belief in smoothness, which is 
in turn reflected in the small diagonal uncertainties. Said an¬ 
other way, if the form factor can be faithfully described with 
a continued fraction with only a few terms, then Re is tightly 
constrained. But does the limited freedom of the continued 
fraction fit inappropriately force a small value of R^? We shall 
also discuss choosing other forms that may drive Re one way 
or the other, especially if we allow undulation in an otherwise 
smoothly falling form factor. We note that our value for Re is 
consistent with the analyses of the Mainz data by Lorenz et al. 
using theoretically motivated analytic forms for Ge Esma. 
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B. GeIGm 

We have extracted Ge assuming that pLpGe/Gm = 1 — 
2 ^/2o 2^ = 8 GeV^, and now wish to investigate conse¬ 

quences of different choices for Ge/Gm, a least to the extent 
of considering other choices for Q^. Fig. |^shows the result¬ 
ing x^/dof for various values of Qq upon fitting the 1422 data 
points with a 4-parameter CF function. The full data set favors 
a value of Qg « 8 GeV^. This is bounded sharply on the low 
side and weakly on the high side. The numbers shown beside 
each point are the values of the extracted radius Re which are 
only slightly influenced by the Ge/Gm ratio used to extract 
Ge from Eq. (j^. Moreover, the radius Re = 0.84 fm is stably 
reproduced for 4 < Ql< 20 GeV^. In particular, for steeper 
R-pGe/Gm slopes, the extracted Re actually decreases slightly. 
This is a result rather different from Bernauer et ai, who ob¬ 
tain a larger radius and {IpGe/Gm = 1 — 2^/(1-4 GeV^) at 
very low Q^. 



FIG. 4. (Color online) Minimum ;t;^/dof values for continued frac¬ 
tion fits of Ge extracted using ^ipGE/G m = 1 — 2^/2o '''i^h different 
values of Qq, when using the Mainz data. Each point shows the value 
of Re extracted for a given Qq. The quantity Re remains stable over 
the wide range 4<2g< 20 GeV^, indicating that Re is not sensitive 
to the Ansatz for Ge/Gm- 

Fig. [5] s hows the world’s polarization transfer data 
l25l I28H33I for iXpGe/Gm on the proton (over a somewhat 
wider range than we have considered in the bulk of this pa¬ 
per). It also shows the fit we have used, and two other fits. In¬ 
cidentally, fitting the form f{Q^) = 1 — 2^/2o 1*^®^ the data 
gives 2 q = 8.02 ±0.05 with j^/dof = 2.3. Although the re¬ 
coil polarization method is the best way we know to determine 
the electric to magnetic form factor ratio, being relatively free 
of two-photon effects, the data points below « 0.8 GeV^ 
disagree with each other more than their quoted uncertainties 
would allow. The variate m,- = [y; —/(x,)]/(T,- has a small mean 
of 0.02 and a large standard deviation of 4.0. From this we 
conclude that a linear fit acceptably represents the average of 
the data points, despite their underestimated uncertainties. 

That the Mainz data also prefer 2o = 8 GeV^, consistent 
with the recoil polarization data, leads to a conundrum. Earlier 
Rosenbluth results gave scaling. The drop in Ge/Gm with 


increasing was a great surprise when announced in 1999 
and published in 2000 11211 . Why the Mainz data, without 
full hard two-photon corrections, agrees with the polarization 
results is a mystery. 



FIG. 5. (Color online) Fits to world polarization transfer data. 
The data are from GEp I (21]HD, GEp II (22||23l, GEp III l25l . 
Zhan et al. 1301 . Ron et al. 1281 . Crawford et al. 1341 . RSS 1351 . 
Paolone et al. (36), and Strauch et al. (371. Some data with larger 
uncertainty limits have been omitted. The solid orange line uses 
IIpGe/Gm = 1 - Q^/Ql with Ql = 8.02 GeV^, the black dashed 
line is the Bernauer et al. (Ml fit up to about 0.3 GeV^ and a hybrid 
produced by other Mainz workers beyond that 1381 . and the green 
dotted line is a fit from Punjabi et al. 1391 . 


C. Epsilon Dependence 



FIG. 6. (Color online) Ge vs. e. The six sets, from top to bottom, 
correspond to the six beam energies of the experiment, 180, 315,450, 
585, 720 and 855 MeV. 


Fig. 0 shows Ge versus e, for the Mainz data, for varying 
Q^. There are 6 sets of points corresponding to the differ¬ 
ent beam energies of the experiment. Since Ge is a function 
of 2^ (and the Ge obtained from data will reflect this if all 
corrections are made), a horizontal line on this plot intersects 
the values of e represented by points at fixed Q^. Likewise, a 
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vertical line on the plot shows different values of at con¬ 
stant e. The large range in e covered allows us to determine 
Ge / Gm, as discussed earlier. It is also possible that some of 
the apparent e dependence can be attributed to mismatches in 
the normalization of the data from different beam energies. 



FIG. 7. (Color online) Dependence of the extracted Ge values on e 
for 20 bins in 2^. The data show little or no dependence on e within 
statistics. (If desired, the actual 2^ values may be inferred from the 
next Figure.) 

Fig. 1^ shows the data set plotted versus e for 20 bins in Q^. 
To ensure a common for each horizontal line in this plot, 
data within a given bin were evolved to a central using 
the 

G£(2Lrage) = GfCgLasured) ■ (18) 

The lines in Fig. |7] look very flat, and the actual slopes of 
the data in Fig.j^are shown in Fig.[^ On average, these slopes 
are zero, but there are small variations. At low Gm barely 
contributes to the cross section, so any e-slope must be an 
indication of a normalization mismatch in data from different 
beam energies. The slopes are never more than a fraction of 
a percent, which can easily be accounted for by systematic 
variations in the data. 


D. Two Photon Contributions 

A possible cause of real e-dependence stems from two- 
photon exchange effects. Although the Mainz data set in- 



Q2 (GeV2) 

FIG. 8. (Color online) e slopes for 14 bins in 2^- 


eludes Coulomb corrections following McKinley and Fesh- 
bach Il40l . which are for the limit of very heavy pointlike pro¬ 
tons, there are further hard two-photon effects occasioned by 
the hadronic structure of the proton ll4TH43l . 

Recent data on the cross section ratio of positron to elec¬ 
tron elastic scattering from the proton verifies the idea that the 
Rosenbluth extraction of the Ge/Gm ratio receives significant 
corrections from two-photon exchange ll44l |45]| . 

A potential further consequence of two photon exchange is 
that in addition to changing the £ slope in the reduced cross 
section Gm{Q^) + [e/x)G\{Q^), there could also be terms 
quadratic or higher in £ ll46l . However, the data show that 
any terms are small. We conclude that two-photon correc¬ 
tions, although they have been demonstrated to exist Elllsl, 
do not induce strong curvature in the Rosenbluth plot or bias 
the data in such a way as to change the radius Re if one does 
not have data over the full range of £. 


E. Fitting with Polynomials 

The continued-fraction fit has only a few parameters and 
may not accurately describe inflections, if there are physical 
inflections, of the measured form factor. In this subsection, 
we experiment with other fits to the full range of the Mainz 
data. We consider five different generic types of fit functions, 
in addition to the CF fit already presented, and will show some 
comparison of the different fits in Fig.|^ See also B9l - [54l . 

First, for reference, we fit the whole data set to a dipole 
form ao(l +2^/«i) Although the x2/(jof is larger than 
acceptable at 2.28, the fit visually is remarkably good, and 
Re = 0.8299 ±0.0002 fm. The dipole is famous for giving a 
small radius when fit to a long data set, although this value is 
intriguingly close to the muonic hydrogen Lamb shift value of 
0.841 fm. 

Second, we have followed the lead of Bemauer et al. a 
and fit to a double dipole f{Q^) = + Q^l‘^ 2 )^^ ± 

(1 — ai)(l ± Q^/a^)^^]. The fit has x^/dof=1.6—the best ap¬ 
parently we can achieve with a smoothly and monotonically 
falling fit function—and Re = 0.859 ± 0.001 fm. 
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Third, we consider polynomial fits. We do not advocate us¬ 
ing polynomial fits beyond the spacelike reflection of the nn 
threshold, since convergence of the fit is not assured beyond 
this point. However, they have been used elsewhere, and we 
would like to comment on the results. Polynomial fits with 
sufficient terms offer flexibility to fit inflections in the data, 
but they inevitably diverge outside any fit region, and accu¬ 
racy at the end points is often poor for global fits. 

Fourth, we consider inverse polynomials ao/{l + 

Fifth, we consider power series expansions in z{Q^), f{z) = 
ao(l +'Lf^^aiz‘), as advanced in this context by Hill and Paz 
ll55l , where 


a/ 4ml + Q'^+ 2mx 


(19) 


the precise behavior of the fit function at the origin signifi¬ 
cantly influences the extracted value of Re- 
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FIG. 9. (Color online) Differences between various fit forms and the 
standard continued fraction fit. All the power-series fit forms show 
undulations when enough parameters are included. The ;f^/dof drops 
accordingly to about 1.37, but at the expense of pathological behavior 
at the origin and above = 1 GeV^. 

Fig. shows a visual comparison of five fits: the dipole fit, 
the double dipole fit, and representatives of the other three fit 
types. The curves correspond to the differences between each 
model tested and the 4-parameter continued fraction (CF) fit 
described earlier. All polynomial fits show multiple oscilla¬ 
tions around the CF value. Moreover, the curves are clearly 
pathological near = 0 and = 1 GeV^, that is to say, just 
outside the region where the fitted data has support. With suf¬ 
ficient parameters, the polynomial, inverse polynomial, and z- 
fits all start to reproduce inflections in the data, and they track 
each other roughly. The large rise at 2^ = 0 is the reason these 
fits give a larger radius than the CF fit. Although in absolute 
terms, the fits differ from each other by less than the point-to- 
point uncertainties on the data, and absolutely less than 0.001, 


FIG. 10. (Color online) Renormalization constants for the 34 Mainz 
normalization sets. Red inverted triangles, black dots, and purple 
upright triangles correspond to Spectrometers A, B, and C, respec¬ 
tively. The numbers are the 6 beam energies in MeV, and the points 
to the left of the corresponding arrow are the sets at that energy. The 
average and standard deviation of these normalization constants are 
1 and 0.015%, respectively. 


F. Systematic Deviations Between Fit and Data 


The typical quoted uncertainty on each point is a few tenths 
of a percent. Since these data represent 34 separately normal¬ 
ized data sets taken with three spectrometers, it is not unrea¬ 
sonable to suppose that some of the apparent undulation could 
be modified or removed by relative renormalization. Since ab¬ 
solute normalizations in each spectrometer are not known to 
better than a percent, there is some freedom to do this on the 
level of at least a few tenths of a percent. 

These sorts of relative renormalizations were made by 
Bemauer et al, with slightly different renormalizations for 
the different fit functions they used. We did a similar pro¬ 
cess using the continued fraction fit. For each normalization 
set we formed the uncertainty-weighted average of the ratio 
Ge{^)/ f{^) for all points in each subset. The data were 
then divided by the inverse of this ratio. These factors are 


shown in Fig. 10 for the various data sets. Arrows indicate the 


beam energy of the points to the left of the arrow. The over¬ 
all renormalization is unity, with a point-to-point variation of 
about 0.15%. This indicates that the original normalizations 
were done well, although they could be be modified a bit when 
using a different fit function. 

Attempting make it easier to see any systematic effects 
within this thicket of renormalization ratios, we combined 
points from different spectrometer settings within the full data 
sample into 14 bins of 2^ with roughly 100 points in each 


bin. Fig. 11 shows the results of this exercise before the 
renormalizations were done (top) and after (bottom). Each 
plot shows the average differences, Ge{Q^) — f{Q^)- There 
are what appear to be statistically significant variations in 
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SDectrometer A 



FIG. 11. (Color online) Average differences {Ge)! — f(Qj) for 14 
bins. Upper: before renormalizing the 34 data sets; lower: after 
renormalizing. 

Ge{Q^) — f{Q^) in the upper plot. In fact, focusing on the 
low-2^ behavior, the data show a trend favoring a larger slope, 
and a bigger radius Re, than our fits suggest. However, it is 
worth noting that these variations are on the order of 0.1%, 
which is commensurate with the point-to-point uncertainties. 
In the lower plot, the systematic variations seen in the upper 
plot are reduced by half, but renormalization cannot account 
for the remaining fluctuations which are on the order of 0.001. 

; separates Ge{Q^) — f{Q^) into plots for each spec¬ 


Fig. 


12 


trometer individually. Here there is a gradual rise and fall— 
albeit on the level of a tenth of a percent—in Ge (Q^) —f{Q^) 
for Spectrometers A and C. For Spectrometer B, which is the 
workhorse at low Q^, there is no such variation. Any devia¬ 
tions in the data from the continued fraction fit should show 
up in all three spectrometers if they are real. Because this is 
not the case, the observed fluctuations likely are not intrinsic 
to Ge- 

We remind ourselves that the 4-parameter continued frac¬ 
tion fit to the full 1422-point Mainz data set is good, and it 
is now somewhat improved by the renormalizations. After 
the renormalizations. Re does not change appreciably, but the 
X^/dof decreases to about 1.4, which by some measure is still 
too large. Consequently, we need to consider increasing the 



FIG. 12. (Color online) Differences, Ge{Q^) — f[Q}), averaged over 
14 2^ ranges for Spectrometers A, B, and C individually, after renor¬ 
malizations of the 34 subsets. 


size of the uncertainty limits. 


V. FINAL RESULTS 


The systematic deviations from the CF fit shown in Fig. 12 
differ considerably from spectrometer to spectrometer, sug- 
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gesting that they are not intrinsic to Ge{Q^) and perhaps 
that the point-to-point systematic uncertainties are underes¬ 
timated. Bernauer et al. themselves have rescaled the un¬ 
certainties per normalization set by factors ranging from 1.07 
to 2.3 Q. Therefore, we repeated this exercise globally for 
Ge{Q^), and found that the uncertainties required rescaling 
by a factor 1.15. Fig. 13 shows the full data-set including our 
renormalizations of data sets and uncertainties. The resulting 
new fit (Fig. 131 has a j^/dof of unity for the full data set. 
All of the modifications we have made to the data have not 
changed Re more than a few parts per thousand. We obtain the 
value Re = 0.8404 fm, from this procedure, with a diagonal 
uncertainty of 0.00007 fm. The value of Re fm remains con¬ 
sistent with the muonic hydrogen Lamb shift measurements. 



FIG. 13. (Color online) Final global fit with the 34 data sets renor¬ 
malized and the point-to-point uncertainties on Ge scaled up by 15%. 


Regarding the size and distribution of the uncertainties, the 
across-the-board increase of the uncertainty limits on Ge by 
15% yields a normal distribution for [GEiQj) — f{Qj)]/(Ji. 
Fig. shows the resulting histogram of this quantity for all 
1422 points in the data set. A Gaussian fit yields a mean of 
zero and a standard deviation of 1 with a good j^/dof, as ex¬ 
pected for Gaussian statistics. 

The statistical uncertainty from the fit is small, and the 
overall uncertainty on Re is dominated by systematics. We 
have estimated the systematic uncertainties by finding the 
spread among a set of extracted radii; 1) Spectrometer B, 
< 0.02 GeV^, fits to terms up to constrained by as¬ 
suming an exponential, Gaussian, or empirical charge distri¬ 
bution: 0.836, 0.849, and 0.859 fm, respectively; 2) < 1.0 

GeV^, fits to the double dipole, continued fraction, and in¬ 
verse polynomial fit forms: 0.830, 0.840, and 0.870 fm, re¬ 
spectively; and 3) a global fit to the 1422 points with each 
of the 34 normalization constants as free parameters, 0.827 
fm (not reported in detail here). The average and standard 
deviation within this set are 0.844 fm and 0.016 fm. We 
take this standard deviation as an estimate of the systematic 
uncertainty on Re from fit-model dependence, and keeping 
the central value from our previous analysis, conclude that 
Re = 0.840 ± O.OOlstat ± 0.016syst- 


c 

3 

O 

O 


-r- 


Entries 

1422 


Mean 

-0.0003825 


RMS 

1.104 


! ndf 

42.27 / 66 


Constant 

51.12 + 1.76 


Mean 

0.004946 + 0.029582 


Sigma 

1.082 ±0,024 




-2 0 2 
[GE(Q2>-f(Q0]/a 


FIG. 14. (Color online) Distribution of [GEiQ]') —/(2?)]/ct/ for all 
data points. Flere the individual uncertainties on each point, C7/, have 
all been rescaled by a factor of 1.15. 


VI. CONCLUSIONS 

The Mainz data set is of extremely high quality— 
expansive, accurate, and self-consistent. We began with an 
analysis of the low-Q^ part of the data set, where a poly¬ 
nomial expansion of the form factors should converge, and 
which should and does yield an accurate result for the pro¬ 
ton radius. We found a proton radius in agreement with the 
muonic hydrogen Lamb shift results and significantly smaller 
than the CODATA value. 

We also analyzed the full data set, assuming that Ge is 
monotonically falling and inflectionless, and used a contin¬ 
ued fraction form to map this. We rescaled the different data 
sets on a level that is smaller than the original normalization 
uncertainties. We inflated the point-to-point systematic un¬ 
certainties by 15%, which is well within reasonable system¬ 
atic uncertainties for such an electron scattering experiment. 
We can then fit all data nicely using only 4 parameters. This 
results in a j^/dof of unity, and with some further consider¬ 
ation of other ways to fit the data, determine a proton radius 
/?£ = 0.840 ±0.016 fm. 

This result is in excellent agreement with the muonic Lamb 
shift results. Of course, if this is true, the proton is less inter¬ 
esting than we had hoped, and beyond-the-standard-model ex¬ 
planations of the proton radius puzzle will not be needed. One 
way to solve the current conundrum using scattering data is to 
have independent confirmation of the shape of Ge from other 
electron or muon scattering measurements, and a number are 
underway or in planning Il56l45^ . In addition, a host of other 
relevant experiments are also underway or under analysis, in¬ 
cluding the completed but not yet published measurements of 
nuclear radii in other muonic atoms ll60l . and new high pre¬ 
cision atomic level splitting experiments that will yield new 
and precise measurements of the proton radius |ImT 4641 . We 
eagerly await all the new measurements that can elucidate the 
proton radius quandary. 
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